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Nuclear  magnetic  resonance  spectrum  of 31 P donors  in  silicon  quantum 
computer 

A.  A.  Larionov,  L.  E.  Fedichkin,  A.  A.  Kokin  and  K.  A.  Valiev 

Institute  of  Physics  and  Technology,  RAS 

34,  Nakhimovsky  pr.,  Moscow,  117218,  Russia 

Abstract.  The  influence  of  an  electric  field  created  by  a gate  potential  of  the  silicon  quantum 
computer  on  the  hyperfine  interaction  constant  (HIC)  is  obtained.  The  errors  due  to  technological 
inaccuracy  of  location  of  donor  atoms  under  a gate  are  evaluated.  The  energy  spectra  of  electron- 
nuclear  spin  system  of  two  interacting  donor  atoms  with  various  values  of  HIC  are  calculated.  The 
presence  of  two  pairs  of  anticrossing  levels  in  the  ground  electronic  state  is  shown.  Parameters  of 
the  structure  at  which  errors  rate  can  be  greatly  minimized  are  found. 


Introduction 

One  of  the  most  perspective  ideas  on  experimental  realization  of  the  quantum  computer  is 
the  possibility  of  creating  artificial  multispin  systems.  Such  an  approach,  not  yet  realized, 
was  offered  by  Kane  [1  ] and  was  studied  further  in  [2—4] . The  special  feature  of  this  variant 
consists  in  the  realization  of  the  individual  addressing  to  separate  spins-qubits.  For  this 
purpose  it  is  supposed  to  use  a silicon  structure  of  the  MOS  type,  where  donor  atoms  of 
stable  isotope  of  phosphorum  31 P replacing  silicon  atoms  in  the  knots  of  a crystalline  lattice 
are  introduced  into  the  thin  layer  of  the  spinfree  silicon  28Si  isotope.  Such  a donor  has 
nuclear  spin  with  7=1/2  and  states  which  are  characterized  by  a high  value  of  an  effective 
Bohr  radius.  In  the  considered  approach  it  is  necessary  to  use  rather  low  temperatures  so 
that  the  electrons  of  donor  atoms  occupied  only  the  lower  spin  state  in  the  magnetic  field, 
that  is  T <SC  2/ub  B/k,  where  mb  = 9.27  x 10_24J/T  is  Bohr  magneton,  B is  induction 
of  the  external  constant  magnetic  field,  k = 1.38  • 10-23  J/K  is  Boltsman  constant.  At 
fields  B > 2T  it  corresponds  to  temperatures  T < 0. 1 K,  which  are  much  lower,  than  the 
temperature  of  freezing  out  the  donors  electronic  state.  Therefore,  the  donors  will  remain 
in  the  neutral  ground  orbital  S-state. 

Each  donor  atom  with  its  nuclear  spin  in  a semiconductor  structure  is  supposed  to  be 
arranged  regularly  with  an  adequate  accuracy  under  “its”  own  control  metal  gate  (gate  A), 
separated  from  the  surface  of  silicon  by  a thin  dielectric.  The  gates  A form  a linear  chain 
of  any  length  with  the  period  /.  The  variation  of  the  electrical  potential  of  gates  J,  located 
between  gates  A,  allows  to  control  a degree  of  overlapping  of  electrons’  wave  functions 
located  on  the  adjacent  donors  a and  b and  the  constant  of  their  exchange  interaction  / as 
well  as  the  constant  of  scalar  interaction  of  their  nuclear  spins  Iab. 

This  is  achieved  by  redistribution  of  electronic  density  between  adjacent  donor  atoms. 
It  is  supposed  that  with  the  help  of  electrical  field  created  by  gates  A it  is  possible  to 
change  distribution  of  an  electronic  density  near  the  nucleus  in  the  main  state  by  adjusting, 
respectively,  individual  resonance  frequency  of  a nuclear  spin  of  each  donor  atom,  which 
is  determined  by  hyperfine  interaction  between  the  atom  and  electronic  spin.  It  allows  to 
realize  quantum  operations  by  selective  influence  of  resonance  radio-frequent  impulses  on 
nuclear  spins  of  certain  donors. 
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Thus,  for  the  development  of  a variant  of  the  semiconductor  quantum  computer  proposed 
in  [1],  it  is  important  to  have  a detail  information  about  properties  of  a donor  atom  energy 
spectrum  hyperfine  structure. 


1.  Dependence  of  the  hyperfine  interaction  constant  on  the  electrical  field 

The  spin  Hamiltonian  of  hyperfine  interaction  is: 

H,s  = A (is)  , (1) 


where  I and  S are  nuclear  and  electron  spins  respectively, 

8"T  Li 

A(V)  = — 1*0(0,  V)|22Mb#nMnP-.  (2) 

3 Air 

/hn  = 5.05  • 1 0“27  J/T  is  nuclear  magneton,  po  = Air  ■ 10-1  T2cm3/J,  = 2.26  is 
Lande’s  factor  of  31P,  hb  = 9.27  • 10~24  J/T  is  Bohr’s  magneton,  4>o(0,  V)  is  the  wave 
function  of  the  electron  near  nucleus.  Taking  into  account  the  existing  electric  field  the 
Schrodinger  equation  for  this  wave  function  should  be  supplemented  by  an  perturbation 
operator  A H = erpip' , z')  (e  — 1.6-  10~19C,  p'2  — x'2  + _y'2)  [2],  where  the  potential  (p  is 
induced  by  the  gate  which  has  the  form  of  round  disk  of  radius  a = l a /2.  By  factorizing  this 
potential  and  taking  as  undisturbed  wave  functions  hydrogen-like  functions,  we  calculate 
matrix  elements  of  transitions  necessary  for  calculation  of  error  corrections  for  the  wave 
function  within  the  framework  of  the  theory  of  perturbations.  Limiting  ourselves  by  the 
second  order  of  the  perturbations  theory  and  conducting  necessary  number  substitutions, 
we  shall  receive  the  following  expression  for  a relative  error  correction  of  the  donor  atom 
hyperfine  interaction  constant  (with  c = 2 a = lOnm,  V is  expressed  in  volts): 


A A(V)  1 

— = 0.55  V - 0.09  V2.  (3) 

A 


Let  us  consider  the  case,  when  gates  A have  the  form  of  endlessly  long  strips  with  the 
breadth  of  2 a — l a located  at  distance  D from  the  conducting  substrate,  from  which  the 
potential  of  a gate  is  accounted.  Conducting  similar  operations,  we  shall  receive  for  a 
relative  error  correction  of  the  donor  atom  HIC  the  following  expression  (with  c — 2 a — 
lOnm,  D/a  — 100): 


AA(V) 

A 


-0.063V2. 


(4) 


The  voltage  in  (4)  is  expressed,  similar  to  the  above,  in  volts. 


2.  Influence  of  the  technological  location  spread  of  donor  atoms  on  the  HIC.  Voltage 
error  on  a gate 


Technological  inaccuracies  of  a donor  atom  location  under  a gate  in  respect  of  the  position 
of  x — 0,  z — c shall  be  designated  by  <5x  and  8z.  We  resrtict  ourselves  here  to  the 
consideration  of  only  the  case  of  a strip  gate.  Factorizing  the  electric  field  strength  and  its 
gradient  along  axes  x and  z to  degrees  8x  and  8z,  we  receive  a resultant  expression  for  a 
relative  error  of  the  HIC: 


2 2c  ) 
2 + c2J 


+ (Sxf 


0.063V2 


2c2  - a2 
(a2  + c2)2 


0.085V 


2c4  - «4  ] 

2c2  ( a 2 + c2)2  j 


(5) 
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Fig.  1.  Levels  of  energy  of  spins  for  a system  of  two  atoms  31 P in  a magnetic  field  ( Aa/J  = 0.3; 
Ah/J  = 0.4). 


The  second  item  in  (5)  can  be  turned  to  zero  with  different  values  of  a,  c and  V.  That  is  why 
by  setting  a required  relative  error,  it  is  possible  to  determine  parameters  of  structure  a and  c 
and  voltage  V at  which  the  required  accuracy  will  be  realized.  Permissible  deviation  of  the 
donor  atom  along  axis  z connected  with  the  technological  inaccuracy  of  location  in  2-3  nm 
corresponds  to  a 1%  preset  relative  error  of  hyperfine  interaction.  Voltage  error  at  a gate 
allowed  for  a correctly  controlled  quantum  computing  process  can  be  determined  from 
resultant  expressions  for  changes  of  the  hyperfine  interaction  constant.  The  bandwidth 
of  resonance  impulse  frequencies  should  be,  at  least,  one  order  less  than  the  resonance 
frequency  of  nuclear  spins  which  is  equal  to  hundreds  of  kHz.  Taking  the  value  of  the 
error  for  the  hyperfine  structure  constant  corresponding  to  the  bandwidth  of  the  resonance 
impulse  <5  (A A)  ~ 104Hz,  we  shall  find  the  value  of  the  error  for  the  voltage  at  the  gate 
being  equal  to  10-4—  10_3V.  Let  us  note,  that  in  case  of  round  disks  gates  because  of 
the  linear  V summand  in  expression  (4),  there  is  a chosen  voltage  at  which  the  allowable 
voltage  error  greatly  increases. 

3.  Energy  spectrum  of  an  electronic-nuclear  spin  system  of  two  interacting  donor 
atoms 

Let  us  consider  two  donor  atoms  being  at  a distance  of  l from  each  other.  This  distance 
should  be  such  that  constant  ./  of  the  effective  exchange  interaction  of  electrons  caused 
by  partial  overlapping  of  their  wave  functions,  had  the  magnitude  sufficient  for  organizing 
double  qubit  operations.  The  spin  Hamiltonian  of  such  a system  looks  as  follows: 

H — 2)j,b B (sa  + + J (sjb)  — ( ' Ja  + h ^ + Aa  ^ IaSa ^ + Ab  ^ IbSb ^ (6) 


where,  Sa,  I a,  Sb , h are  spin  operators  of  both  an  electron  and  a nucleus  for  the  first  and 
second  atoms  respectively,  B is  magnetic  field  (it  is  directed  along  the  z axis),  Aa  and  A/, 
are  constants  of  hyperfine  interaction  depending  on  the  potentials  on  gates,  n\>,  is  Bohr 
magneton,  fi n is  nuclear  magneton,  gN  is  Lande’s  factor  for  ' 1 P.  The  eigen  values  of  a 
Hamiltonian  was  accurately  calculated.  With  the  help  of  this  matrix  the  energy  levels  were 
numerically  established. 
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4.  The  conclusion 

The  proposed  paper  has  a detailed  analysis  of  the  influence  of  the  electric  field  created  by 
gate  A on  the  hyperfine  interaction  constant  for  two  types  of  gates:  in  the  form  of  a disk  and 
in  the  form  of  a strip.  A weaker  dependence  of  the  hyperfine  interaction  constant  on  the 
depth  of  the  donors  position  under  a gate  as  well  as  decreasing  of  its  value  as  the  distance 
between  the  gate  and  the  substrate  was  increasing  were  noted.  The  conducted  calculations 
with  the  two  types  of  gates  make  it  possible  to  conclude  that  there  is  a significant  dependence 
of  the  hyperfine  interaction  constant  on  the  form  of  a gate. 

The  obtained  expressions  enable  to  detennine  permissible  values  of  technological  in- 
accuracies Sx  and  Sz  with  the  set  value  of  the  A constant  error  and  working  values  of  the 
V gate  potential. 

The  paper  contains  calculation  of  the  full  energy  spectrum  of  an  electronic-nucleus  spin 
system  of  two  interacting  donor  atoms  for  different  values  of  HIC.  The  presence  of  two  pairs 
of  anticrossing  levels  in  the  main  electronic  condition  corresponding  to  the  projection  of 
full  electronic  and  nuclear  spins-1  (dotted  line)  and  0 (dashed  line)  is  shown.  For  organizing 
computing  and  measuring  process  both  pairs  anticrossing  levels  can  be  used. 
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